Coupling of phonons and spin waves in triangular antiferromagnet 
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We investigate the influence of the spin-phonon coupling in the triangular antiferromagnet where 
the coupling is of the exchange-striction type. The magnon dispersion is shown to be modified 
significantly at wave vector (2tt, 0) and its symmetry-related points, exhibiting a roton-like minimum 
and an eventual instability in the dispersion. Various correlation functions such as equal-time 
phonon correlation, spin-spin correlation, and local magnetization are calculated in the presence of 
the coupling. 
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I. INTRODUCTION 

A number of recent experimental breakthroughs has 
revived interest in the phenomena of coupling between 
magnetic and electric (dipolar) degrees of freedom in a 
class of materials known as "multiferroics" 1 . Some note- 
worthy observations include the development of dipole 
moments accompanying the helical spin ordering 2,3 , dis- 
placement of magnetic ions at the onset of magnetic order 
in the triangular lattice YMnOs 4 , and adiabatic control 
of dipole moments through applied magnetic fields 5 ' 6 , all 
of which unambiguously point to the strong coupling of 
electric and magnetic degrees of freedom. A number of 
theories has been advanced to establish a microscopic 
understanding of these couplings 7 13 . 

The known mechanisms of the spin-polarization cou- 
pling fall into two categories. One is of the inverse 
Dyzaloshinskii-Moriya (DM) type, whereby the local 
dipole moment, denoted by Uij : couples to the spins Si 



by 



&ij x (Si x Sj). Here the unit vector e^- con- 



nects the centers of the magnetic ions at i and j. The 
microscopic origin of such coupling was investigated in, 
for instance, Refs. 7,11. The behavior of a large class 
of multiferroic materials can be understood on the basis 
of this coupling 1 . The other type arises from exchange- 
striction, wherein the movement of the magnetic ions is 
assumed to directly influence the exchange integral and 
lead to the coupling ~ e^- • (ui — Uj)Si • Sj. The spin- 
lattice coupling in RM^Os (R=rare earth) is believed 
to arise from this mechanism 14 . 

More recently, the dynamical aspect of the mag- 
netoelectric coupling has been investigated both 
experimentally 15 ' 16 and theoretically 17 . The dynamical 
coupling is important because it can arise without the or- 
dering of one or both of the degrees of freedom, and can 
substantially influence the ac dielectric response 15 17 , or 
even result in an exotic new phase with vector chirality 18 . 
The dynamics of the small fluctuations in the ordered 
phase of both the polarization and the spin were exam- 
ined in Ref. 7 for the one-dimensional frustrated chain. 
A corresponding analysis of the coupled fluctuations has 
not yet been tried in the case of the exchange-striction 



mechanism, or for other lattice geometries. 

Triangular geometry offers a potentially fertile ground 
for the interplay of spin and dipolar degrees of freedom 
because of the tendency of spins to form a spiral (120°) 
structure even without further frustrating interactions. 
The ground state is characterized by non-zero (Si • Sj) as 
well as (Si x Sj), which may result in spin-dipole interac- 
tions of both DM and exchange-striction types. Further- 
more, quite recently, it was shown that a spin S = 5/2 
triangular antiferromagnet in RbFe(Mo04)2 19 develops 
spontaneous polarization along the c axis as the spins or- 
der in the ab plane. This and another triangular lattice 
compound, YMnOs, offer promising examples where the 
interplay of spin and dipolar degrees of freedom can be 
revealed in detail. In particular, the dynamical aspect 
of the spin-dipole coupling in the triangular lattice re- 
mains largely unexplored and a theoretical consideration 
of their interplay would be timely. 

In this paper, we examine the coupled dynamics of 
Heisenberg spins and the local dipolar variable (hereafter 
referred to simply as phonons) on the triangular lattice, 
assuming the exchange-striction interaction. In Sec. II, 
the model Hamiltonian is introduced and solved within 
Holstein-Primakoff theory. A number of physical quanti- 
ties, such as the local magnetic moment, phonon correla- 
tion function, and the dynamic spin-spin correlation, are 
derived in Sec. III. Conclusions and the relevance of our 
work to existing experiments can be found in Sec. IV. 



II. SPIN-PHONON MODEL 

The spin-phonon coupled Hamiltonian in the 
exchange-striction picture reads 20 



H- 



(ij) 



Jieji-(uj-Ui)]Si-Sj 



-E 



EL. 

2m 



K 



(1) 



where the antiferromagnetic exchange integral Jij con- 
necting the two nearest-neighbor Heisenberg spins is ex- 
panded to first order in the displacement Ui of each 
atomic site i. eji is the unit vector extending from i 



2 



to j atomic sites. The terms proportional to J\ define 
the spin-phonon coupling. The Heisenberg spin of mag- 
nitude S is represented by Si, and the two-dimensional 
displacement vectors and their canonical conjugate op- 
erators by Ui and pi. We separate the Hamiltonian into 
two parts, H = Hq + Hi, where Hi is the spin-phonon 
interaction term and Ho consists of the Heisenberg and 
phonon Hamiltonians 



(ij) * 

Hi = Ji • (ui — Uj)Si • Sj. 

(ij) 



(2) 



3J S 



V(l-7fe)(l + 2 7fc ). 



(4) 



Here >y k = (1/3) J2 a =i cos l k ' e<*. 
(-1/2, a/3/2), and e 3 = (-1/2, - 
stant is taken to unity. Phonon 
2: directions are also introduced 
well as the phonon energy c^o- 

The spin-phonon Hamiltonian 
second order in the magnon and 
full spin-phonon Hamiltonian to 
in the simple form (k = —k) 22 



with ei = (1, 0), 62 = 
a/3/2). The lattice con- 
operators in the y and 
above as b ky and ft^ as 

Hi can be expanded to 
phonon operators. The 
quadratic order is given 



The classical ground state of the above Hamiltonian was 
worked out in Ref. 20. It was shown that, despite the 
spin-phonon coupling term, the classical spin configura- 
tion is that of the pure Heisenberg model on the trian- 
gular lattice with the neighboring spins at a 120° angle 
with each other. 

Although the spin-phonon interaction does not pro- 
duce any observable effect in the ground state spin con- 
figuration, the excitation spectra of the lattice (phonons) 
and the spins (magnons) will be mixed due to H\. 

The small fluctuations near the ground state can be 
analyzed within the standard Holstein-Primakoff (HP) 
approach after first rotating the spin operators accord- 
ing to their classical spin orientations defined by (Si) = 
(0, sin^, cos fa), where fa is the angle the spin makes 
with the z-axis at site i. All the spins are assumed to lie 
in the 7/z-plane, which also coincides with the plane of the 
lattice itself. After performing the Bogoliubov rotation 
defined by the angle tanh^ = 37^/(2 + 7&) to obtain 
the spin wave spectrum, the Hamiltonian Ho reads 

Ho = Yl £ ka{a k + uo J2 ( b l y b ky + b { z hz) (3) 

k k 

with the spin wave dispersion 



k I 



£ka{ak + uoPkPk + i^k (a\ ~ otj^j (j3 k + /?|) 



(5) 



where 



Xky 
VSXkz 



-;JiS\ 

4 V rnujQ 



e2 \[xt, ' *L 



cos(k • £3) — cos(k • 62), 

2 cos(k • ei) — cos(k • £2) — cos(k • 63). (6) 



Note that only the following linear combination of the 
phonons participate in the interaction with the magnons: 



Xkybky + Xkzbkz 



X k y Xkz 



(7) 



The rotation of the Hamiltonian to the diagonalized basis 
is effected by a series of canonical transformations given 
by ipk = WkXkY^k where ^3/^, the diagonal operators, 
are given by (Aik, A^-, A2k, A^-) T . The respective ma- 
trices are defined as follows: 
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-i sin 
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2A fc 



y/{e k -uj)2+4\l 
v /( £fc -a;) 2 +4A2 



f sinhz^ 



cosh Vk 



a/^Fa 

r fc 

COSh/ii/e = cosh/i 2 fc = ^ 



sinh/ii fe = jjfa, sinh/i 2 /c = 



— &2k 

2k 



(8) 
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FIG. 1: (color online) (a) Dispersion along the (k z ,0) direc- 
tion for Eik (solid) and E 2 k (dashed) for two choices of spin- 
phonon coupling, J\ =0.1 (red and blue lines) and Ji = 0.32 
(pink and violet lines). We have chosen Jo = 3.7 to normal- 
ize the maximum energy value to one. Other parameters are 
S = 1/2, ouo = 0.5, and m = 2. The phonon wave function 
width is chosen 1 / ^muo = 1, equal to the lattice constant. 
The level repulsion is particularly severe at (k z , k y ) = (27r, 0) 
as the strength of the coupling is increased, (b)-(c) Contour 
plots of the low-energy branch E 2 k for (b) J\ =0.1 and (c) 
0.32. Indicated as white dots in (c) are the k points where 
E 2 k reaches zero at the critical spin-phonon coupling. 




FIG. 2: (color online) The new emergent spin configuration 
for the critical spin-phonon coupling value Ji = Ji c corre- 
sponding to the ordering wave vector Q' = (4tv/3 =b 27T, 0), or, 
equivalently, Q' — ±(27r/3, 0). This configuration becomes 
degenerate with those at Q — ±(47r/3, 0) when J\ equals Ji c . 



0.5, and the phonon wave function width mcoo = 1, E 2k 
touches zero at k = (27r,0), (0,27r/\/3), and all their 
sixfold symmetry-related points indicated as white dots 
in Fig. 1 (c). The original zeros of the magnon spectrum 
at ±(47r/3,0) remain intact through nonzero J\. 

As this happens, one has a new spin-ordered pattern 
illustrated in Fig. 2 becoming degenerate with the orig- 
inal 120° ordered phase. This new pattern is obtained 
by rotating spins counterclockwise by 120° along the e\ 
direction and by 60°, also counterclockwise, along the e 2 
direction. The state where the spins are rotated clock- 
wise in both directions will also be degenerate, carry- 
ing the opposite sense of chirality. In terms of ordering 
wave vectors, the new ground states are characterized by 
Q' = ±(2tt/3, 0) instead of Q = ±(4tt/3, 0) as in the 120° 
ordered phase. 



where A k 



= Afccosflfe, T k = A fc sinfl fc , A k = Alfc + A2fc , 



nk,2k 



Sk 

The final form of the Hamiltonian is 

Y£xlwlHkW k X k Y k = dmg(E lkl E 2kl E 2k , E lk ) 
H = Y, (E lk A\ k Ai k + E 2k A\ k A 2 k) • (9) 

k 

The sum is over the entire Brillouin zone. For any 
given k we have E\ k > E 2kl constituting an upper and 
lower branch of the spectra. 

A plot of and E 2k in Fig. 1 shows the change in 
the magnon and the phonon spectrum as Ji is increased. 
The most notable feature in the coupled energy spectrum 
is the appearance of the roton-like minimum at a set 
of /c-points in the Brillouin zone. When Ji equals the 
threshold value, e.g. Ji c = 0.321 for S = 1/2, ujq = 



III. PHYSICAL QUANTITIES 

The local staggered magnetization (uniform magne- 
tization in the rotated basis), (S) = (1/N) J2i(Sf} : is 
modified due to the spin-phonon coupling. The quantum 
correction, defined as the difference of the classical and 
quantum averages S — (£), reads 



■ cosh (f) k (ala k ) - sinh (j) k {a\a\ ) J ,(10) 



which is plotted at T = in Fig. 3 for various spin values 
as the coupling strength is varied. In small powers of 
one obtains the following perturbative expression as the 
quantum correction 



sinh 



2 Wk_ 

2 



^ A| cosh cj) i 



■+E 



A| sinh cj) k 
e k {e k +cj )' 



(11) 
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The first term is the usual quantum fluctuation correc- 
tion, and the latter two account for further corrections 
due to spin-phonon coupling. There is only a tiny change 
in the local magnetization affected by the spin-phonon 
coupling. On the other hand, the upturn found in Fig. 3 
as J\ is driven up to its critical value is probably in- 
dicative of the diverging quantum correction as the new 
ground state is approached. Due to the finite phonon 
mass c^o, the spin-phonon coupling effects are not appar- 
ent until very near the criticality. 
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FIG. 4: (color online) (a) Log plots of the correlation function 
\(m -Uj) — (ui -Uj)o \ in the z and y directions. Using the same 
parameter values as in Fig. 1, the correlation length can be 
extracted as ~ 1.36 and ~ 3.16 lattice constants in the z and 
y directions, respectively, (b) Plot of Gk — 1- Bright regions 
indicate peaks in Gk — 1. 



The equal-time phonon correlation function (u$ • Uj), 
which will be short-ranged (u^ • Uj) = 5ij (since we have 
chosen moo® = 1) in the absence of spin-phonon coupling, 
now reads 



( U( ;.u ) }-(u t .u,)o = ^^*(n- 



■i) 



(G k -l)d 2 k, 



Gk 



■plM+p-k)), 



(12) 



FIG. 3: (color online) Plot of the quantum correction versus 
the coupling strength for various spin values and < J\S < 
JicS. The critical coupling strength depends on S, while the 
product JiS is nearly independent of S. Choices of other 
parameters are the same as in Fig. 1. 
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which is plotted in Fig. 5. The functions appearing in 



at zero temperature. A log plot of \(ui • Uj) — (u* • Uj)o| 
is given in Fig. 4, showing an exponential decay with 
a correlation length which depends on the parameters. 
The momentum space correlation Gk shows pronounced 
peaks around k = (0,27r/\/3) and other symmetry- 
related points as shown in Fig. 4 (b). These are the same 
points where E 2 k shows pronounced minima for large J\. 
Detection of such peaks in the phonon structure factor 
Gk will be instrumental in identifying the spin-phonon 
coupling in a triangular antiferromagnet. 

The spin-spin correlation function, (Sj~(t)S~ (0)), can 
be an effective probe of the spin-phonon coupling. Us- 
ing the straightforward algebra we have calculated the 
absorption spectra I(k,ui) as the imaginary part of the 
Fourier transform of the spin-spin correlation function, 



(13) 



Eq. (13) are defined as = k± (47r/3,0), and 

A\k = sinh [iik (cos Ok + cosh Vk ) — cosh \i\k sin Ok sinh Vk 
A2 k — si nn M2/c ( cos Ok — cosh Vk ) — cosh ji2k sin Ok sinh Vk 
Bik = cosh jiik (cos Ok + cosh Vk ) — sinh \i\k sin Ok sinh Vk 
E>2k = cosh (cosh Vk — cos Ok ) + sinh ji2k sin Ok sinh z/& . 

(14) 



ttS 



((fli fca - A lka )5(w - E lka ) + (B 2ka - A 2k J5(tu - E 2k S) 



-a=l,2 



+ 2e"* fc ((^ifc + B lk )5{u - E lk ) + (A 2k + B 2k )d(co - E 2k ) 
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FIG. 5: Plots of the spectral function I(k,u) along the /in- 
direction (k y = 0) for (a) J\ — 0, (b) J\ — 0.2, and (c) J\ — 
0.32. Emergence of a new low-energy feature at k z = 2tt/3 
for Ji close to the critical value J\ c — 0.321 is apparent in (c). 
(d) Plot of I(k z ,ky,cu = 0.1) clearly indicates new intensity 
peaks at (27r/3, 0) and other symmetry points (elongated and 
shaded) while the bright patterns at (47r/3, 0) and etc. are 
due to the original spin waves. 

The flattening and the eventual collapse of the excita- 
tion band found earlier now manifests itself as intensity 
patterns at (k z ,k y ) = (27r/3,0) and its rotational coun- 
terparts as can be seen in Fig. 5 (d). 

IV. DISCUSSION 

In summary, we have considered the magnon-phonon 
coupling in the exchange-striction coupled triangular lat- 



tice antiferromagnet for Heisenberg spins in the Holstein- 
PrimakofT approach. The dynamics of the lattice and the 
spins are coupled to produce interesting modifications in 
the excitation spectra, in particular (i) the significant 
lowering of the magnon excitation energy at wave vector 
transfer ±27r as indicated in Fig. 1, and (ii) the concor- 
dant variation in the phonon structure factor as shown in 
Fig. 4 (b). Detection of an additional low-energy spectra 
in the spin spectral function I(k,uu) and in the phonon 
structure factor Gk through neutron scattering experi- 
ments will be a clear hint of the strong spin-phonon cou- 
pling. 

Naively speaking, integrating out the phonon coor- 
dinate from Eq. (1) would generate the effective inter- 
action, ~ ~Z)i(Z)j=NN of itijSi ' S j) 2 i which embodies 
the ferromagnetic biquadratic exchange, —{Si • Sj) 2 , and 
some three-body interactions as well. A quantum spin 
model involving quadratic and biquadratic exchanges 
were considered extensively 23 following the discovery of 
the liquid-like ground state in the triangular antiferro- 
magnet NiGa2S4 24 . The ground state revealed correla- 
tions, dynamic on the scale of ~ 1 ns, of a period-six spin 
orientation (60° angles between nearest neighbors), quite 
unlike the period-three orientation (120° angles between 
nearest neighbors) expected in the triangular lattice. It 
is not easy to reproduce such a spin structure in the mean 
field solution of the spin models considered in Refs. 23. 
In fact, the spin-spin correlation observed in NiGa2S4 is 
almost exactly the one shown in Fig. 2. To correctly ac- 
count for the observed periodicity of spins in NiGa2S4, 
one would have to consider an extended-neighbor inter- 
action as in Ref. 25, or some dynamical consequence of 
spin-phonon coupling as in the present paper. 
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